clear;
clc;

S_11_mag = 0.3252;
S_11_ang = -18.15;
S_12_mag = 0.1789;
S_12_ang = 69.01;
S_21_mag = 2.018;
S_21_ang = 72.76; 
S_22_mag = 0.6569;
S_22_ang = -25.83;

det_mag_1 = S_11_mag*S_22_mag;
det_ang_1 = S_11_ang+S_22_ang;
det_mag_2 = S_12_mag*S_21_mag;
det_ang_2 = S_12_ang+S_21_ang;

[det_mag,det_ang] = addPolars(det_mag_1,det_ang_1,-det_mag_2,det_ang_2)

precalc1_mag = det_mag*S_22_mag;
precalc1_ang = det_ang+(-S_22_ang);

[C_1_mag,C_1_ang] = addPolars(S_11_mag,S_11_ang,-precalc1_mag,precalc1_ang)

precalc2_mag = det_mag*S_11_mag;
precalc2_ang = det_ang+(-S_11_ang);

[C_2_mag,C_2_ang] = addPolars(S_22_mag,S_22_ang,-precalc2_mag,precalc2_ang)

B_1 = 1+(S_11_mag)^2-(S_22_mag)^2-(det_mag)^2
B_2 = 1-(S_11_mag)^2+(S_22_mag)^2-(det_mag)^2

GammaS_mag = (B_1-sqrt(B_1^2-4*(C_1_mag)^2))/(2*C_1_mag)
GammaS_ang = -C_1_ang

GammaL_mag = (B_2-sqrt(B_2^2-4*(C_2_mag)^2))/(2*C_2_mag)
GammaL_ang = -C_2_ang